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THEORETICAL LIFT AND DRAG OF THIN TRIANGULAR WINGS AT SUPERSONIC SPEEDS 

By Clinton E. Brown 


SUMMARY 

A method is derived jor calculating the lift and the drag due 
to lift oj point-forward triangular wings and a restricted series 
of sweptback wings at supersonic speeds. The elementary or 
“ supersonic source” solution of the linearized equation of mo- 
tion is used to find the potential function of a line of doublets. 
The flow about the triangular flat plate is then obtained by a 
surface distribution of these doublet lines. The lift-curve slope 
of triangular wings is found to be a function of the ratio of the 
tangent of the apex angle to the tangent of the Mach angle. As 
the apex angle approaches and becomes greater than the Mach 
angle, the lift coefficient of the triangular wing becomes equal to 
that of a two-dimensional supersonic airfoil at the same Mach 
number. 

The drag coefficient due to lift of triangular wings with, 
leading edges wellhehind the Mach cone is shown to be close to 
that of elliptically loaded wings of the same aspect ratio in 
subsonic flight. The resultant force on wings with leading 
edges outside the Mach cone, however, is shown to act normal 
to the surface and thus an induced drag equal to the lift times 
the angle of attack is obtained. 

INTRODUCTION 

In reference 1, Jones calculated the lift of thin point- 
forward triangular wings for the cases in which the apex angle 
of the wing was very small. It was pointed out that the 
results obtained should be applicable in both supersonic and 
subsonic flight, the criterion for the case of supersonic 
flight being that the apex angle be small as compared with 
the Mach angle of the flow. The present paper, making 
use of less restricted theory, extends Jones’ work to the case 
of triangular wings having large apex angles and traveling 
at supersonic speeds. A recent paper was published by 
H. J. Stewart (reference 2) in which the lift of triangular 
wings has been computed, but the method used appears to 
be entirely different. 

In the present theory, the linearized equation of motion 
was used and the results must therefore be restricted to 
small angles of attack and moderate supersonic Mach 
numbers. The solution which has been found should hold 
good for large values of the apex angle up to and coincident 
with the Mach angle. Jones (reference 3) and Puckett 
(reference 4) have found solutions for the drag of triangular 
wings of small thickness at zero angle of attack. The 
solutions are applicable to wings having the leading edges 
either in or out of the Mach cone springing from the apex 
of the wing. Puckett has pointed out that, for the case 
where the leading edge is ahead of the Mach cone, these 


solutions can also be used to calculate the lift; thus, with the 
present solution, the lift for the whole range of apex angles 
at supersonic speeds may be obtained. The pressure dis- 
tributions and lift-curve slopes obtained in the present 
paper can be used to obtain the lift and drag characteristics 
of a limited series of sweptback wings. The drag due to 
lift of the triangular wing has been calculated and a suction 
force has been found to exist on the leading edge. In order 
to use the suction force, however, it appears necessary to 
provide an airfoil section with a rounded leading edge. The 
author is indebted to Mr. Arthur Kantrowitz of the Langley 
Memorial Aeronautical Laboratory for suggesting the method 
used to calculate the induced drag. 

SYMBOLS 

a angle of attack 

A aspect ratio (j^ 

b maximum span of wing 

j3=VM 2 -l 


c 

tangent of apex angle 

c L 

lift coefficient 

c Dt 

drag coefficient due to lift (ffg) 

Di 

drag force due to lift 

E 

source strength 

€ 

apex angle of wing measured from flight direction 

/(*) 

doublet-distribution function 

F 

suction force on wing leading edge 

_£l 


Xi 


i=f—l 

I 

strength of line doublet 






length of wing or root chord 

L 

lift force 

x= f c 

— — — t anh -1 V 1—0 V d<x 

J-c 

y<7 2 -aVl-0V 

M 

Mach number 


Mach angle fsin -1 
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Ap lifting pressure 

g dynamic pressure ^pF 2 ^ 

R distance along leading edge from wing apex 
p density of the fluid 

S wing area 

s„ distance normal to leading edge 

u velocity increment in ^-direction 

U velocity increment normal toJeading edge 

v velocity increment in ^/-direction (j^ 

V flight velocity 

V) velocity in 2 -direction 

W resultant velocity in 2 -direction created by the 
doublet distribution 

x, y, z coordinates of an arbitrary field point 
x u y u Z\ coordinates of a source or doublet 
<£__ disturbance-potential function 

cf > o potential of a supersonic source 

4 > i potential of a line of supersonic sources 

potential of a line of supersonic doublets 
Subscripts: 

ft normal to leading edge 

A triangular-wing condition 

<» infinite-span or two-dimensional wing condition 

THEORY FOR LIFTING TRIANGLE 

The linearized equation of motion of a nonviscous com- 
pressible fluid may be written 

d^dz 2 dx* {1) 

where is the potential function assumed to represent the 
effect of a small disturbance set up by the body being con- 
sidered. The body in this case is a triangular flat plate 
having its vertex at the center of the coordinate system and 
lying in the xy-plane (fig. 1). The problem is to find a 
solution of equation (1) that will satisfy the known boundary 
conditions which are: (1) that the flow be quiescent ahead 
of the Mach cone and (2) that the flow at the surface of the 
plate be tangent to that surface. Because of the linear 
character of the differential equation (1), more general 
potentials can be built up from simple well-known solutions 
such as the one for a single source. 


Mach cone. 



Fiquee 1.— Coordinate system. 

denominator of the integrand is zero. Physically interpreted, 
the range of integration is from the origin to the last 
source point which can influence the field point. Perform- 
ing the integration yields 

x-poy-pTZ , , x-pay-pTZ 

(1— 0V— /JV)®* V(1 - 1 6 V— fiV) 


If two such source lines of opposite strength are brought 
together from the 2 -direction at the a^-plane while the 
product of source strength and the angle between them is 
kept constant, the potential of a line of doublets in the 
xy-plane at an angle tan - V from the x-axis is obtained. Thus 
differentiating with respect to r and setting r=0 gives 


-Iz? 

'(1-0V) 8 ' 2 


(pZq-ctnlr^) 


where 


where The potential of a line of sources with 

strength proportional to x can be found as follows: 

_ f x ' Exi <foi , 0% 

V1 Jo V(x-x 1 )‘ 2 -/3 , (y-vx0 1 -/S 2 (2-rx 1 ) 2 

where a=—i t=—> and x' is the value of x x for which the 

*i *i 


4 V(l-|3 V^-W+z 2 )] 

and I is_ the doublet strength. Differentiating the potential 
function with respect to z gives the vertical velocity w: 


;(j^-ctnh-‘f) 


'(1— /SV) 8/, \|- 2 — 1 V 

. nr., j3 2 (x-jS 2 gy) Vx’-W+s 2 ) 
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It will be noticed that the line doublet creates a conical 
field as the velocity is only a function of z[x and yjx. Since 
the triangular flat plate is a conical body which creates a 
conical field, an attempt will be made to build up the flow 
about the lifting triangle by a suitable distribution of line 
doublets inasmuch as the addition of two or more conical 
fields having the same vertex always creates another conical 
field. The distribution of line doublets must satisfy the 
boundary conditions at the body surface which may be 
written: 

W = Va (7) 

*-o 

where TT is the resultant vertical velocity of the line-doublet 
distribution. If the distribution of line doublets is /(<r) 

(~o<l<o) ( 8 ) 


and substituting equation (12) in equation (13) and inte- 
grating gives 


L=pVIc r VC* 

(14) 

n L 2WC 
L ' L ~qS~ V 

(15) 


The value of the constant I must be obtained by solving 
equation (8). 

The value of the normal velocity at the plate and hence 
the angle of attack may be found by integration of equation 
(8) and letting 2=0. The integration is involved and the 
method of integration is given in appendix A. The resulting 
expression for IT is obtained as 

(16) 


in which tan -1 C=«, the angle of the leading edge. The 
distribution function f(a) can be found in a rather simple 
way by analogy with the solution for incompressible two- 
dimensional flow about a flat plate. Differentiating equation 
(8) with respect to yjx and setting 2=0 gives 



which is of the same form as the integral equation obtained 
when the incompressible flow normal to a two-dimensional 
flat plate is constructed by a doublet distribution. The 
expression /(<r) for the incompressible case would be 

= ( 10 ) 


From this equation the value of I may be calculated. If 
the numerator and denominator of the integrand are multi- 
plied by jS the resulting integral can be seen to be dependent 
upon only the quantity PC or tan e/tan ju. The value of 
the integral may be obtained easily by making the substitu- 
tion p 2 C i —0 2 o 2 =n 2 and plotting the resultant expression. 
This procedure has been followed for values of PC between 
0 and 1 and the result is given in figure 2. The value of I 
is found to be 



(17) 


where \ is the integral term of equation (16). The lift- 
curve slope is now 


C L _2w 2 C 

a ir+X 


(18) 


That this expression is a solution of equation (8) must be 
verified by substitution in equation (8), inasmuch as equation 
(9) is a divergent integral. This proof is carried out in ap- 
pendix A. The value of the velocity in the ar-direction u can 
now be obtained 


u S-/ M w d ' r <>» 

The integration indicated in equation (11) is presented in 
appendix B. The expression obtained from equation (11) on 
the lifting surface (s=0) gives 


u= ± 


rC 2 ! 




( 12 ) 


the sign of the expression being opposite for the two sides of 
the plate. The result in equation (12) shows that the shape 
of the pressure distribution is independent of aspect ratio. 
The lift for the isosceles triangle with root’chord c r is 


L= 




dir 


As PC approaches zero, X also approaches zero and the lift- 
curve slope from equation (18) approaches Jones’ value. 
(See reference 1.) Equation (18) shows that the lift-curve 
slope is a function of only the apex angle and the parameter 
tan e/tan p. It is interesting to note that mathematically 
there is a finite lift-curve slope at the Mach number 1.0. 



( 13 ) 


tan « 

Figure 2.— Variation of X with m 
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The theory is not valid, however, near M= 1 because of the 
original assumptions used in obtaining equation (1). The 
lift-curve slopes of two triangular wings are plotted in figure 
3 against Mach number. 



Fig ore 8.— Variation with Maoh number of Ceja for two triangular wlng3. 


Rewriting equation (18) in a convenient form gives 


,tan_€ 

/Cx\ ta n f i 

\ a ) a (x+X) -yjM 2 — 1 


(19) 


Now, according to Ackeret’s result (reference 5) the lift- 
curve slope of a two-dimensional flat plate is 



A single curve for all Mach numbers can therefore be plot- 
ted if the ratio tan e/tan p is used for the abscissa and 

is r used as the ordinate. This curve is shown in 

figure 4. 

It can be seen that as the apex angle approaches the 
Mach angle the triangular wing provides the same lift co- 
efficient as a two-dimensional wing at the same Mach number. 

The case of the triangular wing having the leading edge 
ahead of the Mach cone from the apex has been treated in 
reference 4. It was found that the lift coefficient obtained 
is the same as that of a two-dimensional airfoil flying at the 
same Mach number. The curve shown in figure 4 therefore 

becomes flat at values of ; — 6 >1. A typical pressure dis- 

tr&n fX 

tribution over a wing having ^-^>1 is shown in figure 5. 
DRAG DUE TO LIFT 

The tliin-airfoil theory used herein gives the result that 
the resultant force is directed normal to the plate, a result 




quite like that obtained from the tliin-airfoil theory at sub- 
sonic speeds. In the solution for subsonic speeds, however, 
a simple extension described in reference 0 permits calculat- 
ing the force due to the suction on the leading edge. It is 
reasonable to suppose that the same method is feasiblo for 
the triangular wing, as the pressure distribution in the neigh- 
borhood of the leading edge is identical in the limit with 
that for^ two-dimensional flat plate in subsonic flow; that 
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is, the velocities normal to the leading edge are in the form 

( 21 ) 


U =-J= 


where s n is the distance measured normal to the edge and 
■Q is a constant. According to reference 6 the force normal 
to the edge in the direction of the velocity is 


F=pttQ 1 


( 22 ) 


Equations (21) and (22) are, however, based on incompressible- 
flow relations and must be corrected for compressibility. 
A simple extension is found if the well-known concepts of 
the Prandtl-Glauert rule are applied. For the two-dimen- 
sional case it is found that the effect of compressibility on a 
flow having a given vorticity distribution in a plane is to 
r educe th e velocities normal to the surface by the factor 
Vl —Al n where M' is the Mach number of the flow. There- 
fore, if the strength of each vortex element is increased by a 

factor the resultant velocities normal to the surface 

will again be equal t-o those of the incompressible flow. In 
this case, however, the tangential velocities at the surface 
and therefore the forces on the surface are increased by the 

ratio conce Pt i* 3 vvell known in thin-airfoil 

theory where it gives the result that the lift-curve slopes of 
thin airfoils are increased by As the total resist- 

ance for Ii/<1 is still zero, this result indicates that the 
leading-edge suction force has been increased over that of 

the incompressible flow by the ratio ' 

It appears that the incompressible equations governing 
the leading-edge suction force (equations (21) and (22)) must 
be corrected as follows: 

The leading-edge suction force on a two-dimensional plate 
will be 


E n = 


ptG 2 

yi -M n 


(23) 


when the vorticity distribution y at the leading edge or the 
tangential velocities are given by the following relations: 


y = 

J «-»0 


20 

yi=3z* vs 


or 


iwo Vl— M' 2 Vs* 


(24) 


The value of the velocity in the y-direction on the triangular 
wing has been calculated to be 

i A 


v= 


V--®' 


Combining this expression with equation (12) gives 


V^fJ 

= ± V«»V1-M /2 


(25) 


where R=x^l-\-C 1 and M‘ is the component of the hight 
Mach number normal to the leading edge. The normal 
force on a small element dR of leading edge from equation 
(23) becomes 


dF A _ ptCRP^T^M* 
~SR 2 


(26) 


and the force on one edge of the isosceles triangle with root 
chord c T is: 


p-PPCc* (1 + C 2 ) Vl-AP 2 
4 


(27) 


and the force in the flight direction from two edges becomes 

F— 2F n sme 

P^VI+^c^Vl-Af^ 


(28) 


Substituting from equation (14) 


F= 


L'-y/l + U* Vl-AP 2 

La / yi+t? Vl 

1 + - 
T 


(29) 


where 6 is the maximum span of the triangular wing. The 
induced drag or, more exactly, the drag resulting from the 
lift may be written 


D ( —La—F 



VT+g Vl-M^ 


1 +- 

IT 


(30) 


Writing the identity 


yi+^ 


equation (30) becomes 


n_W 

Ut ~ 2 


/ yi— gg \ 

r ) 


(31) 


(32) 
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It will be noticed that this result is identical with Jones’ 
result (reference 1) in the limiting case of 0=0. The in- 
duced drag coefficient is found to be 

[ 2 (i +£)- Vl 3 ^] (33) 

where A is the aspect ratio. Equation (33) indicates that 
the triangular wings can obtain a considerable suction force 
at the leading edge and that the drag coefficient due to lift 
of slender wings is very close to that obtained from elliptical 
wings of the same aspect ratio at subsonic speeds. It should 
be pointed out, however, that as soon as the wing leading edge 
passes through the Mach cone, the possibility of obtaining 
a leading-edge suction is gone and the resultant force must 
become normal to the plate surface. This transition corre- 
sponds quite similarly to the case of a two-dimensional air- 
foil passing through the speed of sound. 

DISCUSSION AND CONCLUSIONS 

The lift at supersonic speeds of triangular wings having 
straight trailing edges has been shown to approach the lift 
of a two-dimensional airfoil as the leading edge approaches 
the Mach cone springing from the apex of the triangle. For 
the case where the triangular wing lies behind the Mach 
cone, a suction has been found to exist on the leading edge. 
In order to utilize this suction force in practice it would 
appear necessary, as in subsonic flow, to provide an airfoil 


section with a rounded leading edge. Triangular wings 
should be capable of higher lift-drag ratios than unswept 
wings at supersonic speeds when operating with their leading 
edges not too far behind the Mach cone; the improvement 
should be due to both reduced wave drag and reduced 
induced drag. 

The lift and drag of a series of limited sweptback wings 
may also be calculated with the method developed. It will 
be noted that the pressure distribution over the triangular 
wing cannot be changed if the trailing edge is cut off from the 
tip to the center line along an angle always greater than 
the Mach angle. This fact arises from the nature of the 
supersonic flow in which disturbances cannot propagate any 
farther forward than the Mach cone from the origin of the 
disturbance. , The aforementioned procedure produces there- 
fore a series of tapered sweptback wings haying pointed tips. 
A new series can also be constructed by cutting off the tips 
along lines having angles greater than the Mach angle. 
In each case the pressures over the remaining portions of 
the wing will be the same as though the cutbacks had not 
been made. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Ya., November 29, 1946. 



APPENDIX A 

CALCULATION OF VERTICAL VELOCITY INCREMENT 
The value of IF is, from equations (6) and (8), 

W—I f C ,/£ 7 2 — /T 2 f ft 1 da r , L_1 .o 2 2 ( g ~~j3 2 gy)V^ ^{y 1 ~f~ z 2 ) | 

Integrating by parts gives 

w - z Zkfl@?S-' v=hr* G^i * 


(Al) 


da (A2) 


The integration of the term under the indefinite integral can be performed by parts to obtain the result 




Vz*— Pilf+z 1 ) 

Fax{x>- 20Y-^) + /?W+ z 2 )! 

1 * tun- 1 C **-0*3*)«r-zy 

VI— <SV 



_ a 2 (x t —^z s )—2xya+y i + £ _ 

pz zt/x 2 -P % 2 +z 2 ) 


Substituting in equation (A2) and rearranging terms gives: 


r c dw 


, /VZ*-| 3 2 (j/ 2 +Z*) 

r " 

f2[r i -(8 2 (2/ 2 +s ! )](7a:-^z*-j5V-2/5 I 2 2 ) 

J-cVl-^VVC*-^ 

^ClXLu S, £.2 

1 (*«-/JW) J 

-c^C 2 -o* 

t^fz 2 — 0V) — 2^0-+^*+ z 2 ] 


da (A3) 

It will be noticed in the preceding operation that all the terms containing (he singularity of the form l/z cancel. If 
3 is made to approach zero and terms are collected, equation (A3) becomes 


v-tjlc f)XcVP=7 fi+Kf) 

Completing the integration of the second terms gives 

^(l) 


da 


(A4) 




Differentiating expression (A5) with respect to y/x and performing the integration gives 


da+I^I 


(A5) 


dTT 




*fnt 


»(f) V-KfJ V-Kf) 


=o 


the vertical velocity is therefore constant over the plate surface and the expression f(a) = -yjC 1 — <? is truly a solution 
of equation (8). It is possible that this solution is not unique; however, other solutions would undoubtedly lead to 
physically impossible conditions. 


103 



APPENDIX B 

CALCULATION OF AXIAL VELOCITY INCREMENT 


The value of the 2 -component of velocity may be written from equation (11) 
This integral can be broken up into two separate integrals, as follows: 


(Bl> 


u- 


~~ 7 








2xya 


V-fl’MzM V) 




da 


( 


(T 8 - 


2 xya . g+£ V 

z 8 — |3 i 2 !l " + 'ar i — /3 2 2V 


(B2> 


The first integral appearing in equation (B2) can be integrated by parts to give, for the first complete term, the integral 


Izxy^j g=g(g+gj co 


(zMV 2 ) 8 


J-c -JC 3 —c 


da 


V X 2 -^3 2 + X 2 -|S S Z*/ 

Evaluation of this integral, which may be found in reference 7, equation (228), gives finally 

■ 


(B3> 


: 2 


f 


1 

xy—izij^—fi 1 (y i +z i ) 

" ” xy+izjx 3 — p i (y 3 +z i ) 

yjc*- 

xy-iz- % /x*—p 2 (y i + a?) 

2 pey+isV**— /3 2 (y 2 +z*) “ 


(B4> 


the sign of the expression being opposite for the two sides of the plate. The second integral term of equation (B2) can 
be integrated by breaking the integrand up into four partial fractions, as follows: 




-(a 8 — p 2 z 3 ) !l da 


4W [s’-W+a*)] 3 [<r 

+ 


xy-iz-yj z 4 — ^(V+g 8 ) 

x *_ j3 2 2 * 


^ . ' - -1 p 


(a^-jSV) 2 

xy-iz-ffi—0 3 (y 3 + g^ 3 


xf—fPz 3 

{tf-Fz^yda (z 2 -^) 2 ^ 


(B5) 


The expressions (B5) may now be integrated (reference 7, equation (207)) giving the expression for the complete 
second term of equation (B2) 


±4 


^l c ,. ^—iz^x 3 -PHy i +^) ~J '^l c ,_ pzy+igVz 8 — g^+z 8 )" 


ixJgVz 8 -^(y !i + g) 


V 


Tl 


‘-PTy 
2(z s — jSV) 

I— f / I 1_ 

z 8 — jS 2 ^ 2 

Combining equations (B4) and (BQ) and setting 2=0 yields for ttjon the surface 


s*— 0 V 

zy— feVz 8 — <3 2 (y 8 +g t y 


zy-ftgVz 8 — g(F+^- 


l/i2 ~xy—iz-4^— I rz Vxy-\-izil7?— |8®'(s/*-7-a*)”P 

.V 6 L z^-jsv ; J v 6 L z 4 -^ _ 


«= ±- 


*JC 8 


'CM £ 


(B6) 


(B7) 


REFERENCES 

1. Jones, R. T.: Properties of Low-Aspect-Ratio Pointed Wings at 

Speeds below and above the Speed of_ Sound. NACA Rep. No. 
835, 1946. 

2. Stewart, H. J.: The Lift of a Delta Wing at Supersonic Speeds. 

Quarterly Appl. Math., vol. IV, no. 3, Oct. 1946, pp. 246-254. 

3. Jones, Robert T.: Thin Oblique Airfoils at Supersonic Speed. 

NACA Rep. No. 851, 1946. 

4. Puckett, Allen E.: Supersonic Wave Drag of Thin Airfoils. Jour. 

Aero. Sci., vol. 13, no. 9, Sept. 1946, pp. 476-484. 

104 


5. Ackeret, J.: Air Forces on Airfoils Moving Faster than Sound. 

NACA TM No. 317, 1925. 

6. Von IT A rm rf n, Th., and Burgers, J. M.: General Aerodynamic 

Theory — Perfect Fluids. Theory of Airplane Wings of Infinite 
Span. Vol. II of Aerodynamic Theory, div. E, ch. II, see. 
10, W. F. Durand, ed., Julius Springer (Berlin), 1936, pp. 51, 62. 

7. Peirce, B. 0.: A Short Table of Integrals. Third rev. ed., Ginn 

and Co., 1929. 



